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The  role  of  the  fluctuating  wall  shear  stress  on  the  noise  gener¬ 
ated  by  incompressible  boundary-layer  transition  was  recently 
analyzed  [G.  C.  Lauchle,  J.  Acoust.  Soc.  Am.  67,  158-168  (1980)]. 
The  solutions  of  that  analysis  suggest  that  the  shear  stress 
mechanism  generates  hydrodynamic  sound  weakly  in  comparison  say 
to  the  sound  generated  by  a  fully  developed  turbulent  boundary 
layer  flow.  Because  boundary-layer  transition  is  a  very  unstable 
and  non-steady  flow,  we  might  expect  that  other  acoustic  mechanisms 
exist.  In  this  paper,  we  examine  one  other;  namely,  that  of  the 
fluctuating  boundary  layer  displacement  thickness  which  is  expected 
to  occur  as  turbulence  bursts  are  created  and  convected  over  the 
surface.  The  analysis  shows  that  a  weak  monopole  sound  source 
occurs,  and  at  low  subsonic  flow  velocities,  this  source  generates 
significant  acoustic  energy.  Comparisons  of  the  theoretical 
predictions  with  available  experimental  data  support  this 
conclusion . 
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lightening  discussion  regarding  the  possible  role  of  fluctuating 
displacement  thickness  on  boundary-layer  noise.  This  work  was 
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INTRODUCTION 

This  author  has  recently  analyzed  the  role  of  the  fluctuating  wall 
shear  stress  on  the  radiated  noise  due  to  incompressible  boundary-layer 
transition  [1].  In  the  development  of  that  solution,  several  types  of 
sources  were  identified  under  the  assumption  of  the  boundary  layer  flow 
being  confined  to  one  side  of  an  infinite  rigid  planar  surface.  Quadru- 
poles  and  octupoles  resulted  from  the  fluctuating  Reynolds  stresses  and 
their  images  while  a  longitudinal  dipole  resulted  from  the  shear  stress 
fluctuations.  Under  the  assumption  of  low  Mach  number  flow,  the  dipole 
was  considered  to  dominate  and  was  hence  analyzed  in  detail.  Upon  com¬ 
paring  predictions  using  the  shear  stress  model  for  the  noise  with  experi¬ 
mental  data,  rather  poor  agreement  was  achieved.  Because  the  computed 
radiation  efficiency  was  also  found  to  be  low,  the  discrepancy  between 
theory  and  experiment  was  not  unexpected.  That  is,  the  experiment  was 
performed  with  a  buoyant  body  that  supports  a  substantial  area  of  fully- 
developed  turbulent  flow  which,  based  on  the  magnitude  of  relative 
radiation  efficiencies,  would  be  expected  to  dominate  the  noise  spectrum. 

Although  controversy  exists  as  to  whether  the  transition  zone  is  an 
important  source  of  hydrodynamic  noise,  it  appears  that  further  analysis 
and  experiment  are  required  in  order  to  resolve  the  issue.  In  the 
tutorial  paper  by  Haddle  and  Skudrzyk  [2],  p.  149,  the  authors  state, 
"...the  region  where  turbulence  starts  is  particularly  unstable  and 
intermittent,  and  must  be  expected  to  generate  more  sound  than  the  rest 
of  the  unit."  From  experiments  this  author  performed  in  The  Garfield 
Thomas  Water  Tunnel  [3],  this  statement  was  partially  substantiated; 
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however,  the  theory  of  Ref.  [1]  does  not  predict  any  excessive  noise  over 
that  of  a  turbulent  boundary  layer  alone.  Realizing  that  boundary-layer 
transition  is  a  very  complicated  type  of  flow,  the  use  of  only  one  mecha¬ 
nism  to  describe  its  overall  noise  production  is  probably  not  justified. 

The  shear  stress  mechanism  certainly  contributes,  but  it  may  not  be  the 
dominant  one. 

In  this  paper  we  will  consider  another  type  of  mechanism;  one  that 
gives  rise  to  monopole  sound  radiation.  Its  existence  is  established 
through  a  qualitative  look  into  the  physics  of  intermittent  boundary- 
layer  flow.  It  is  well  known  from  elementary  boundary  layer  theory  that 
the  displacement  thickness  represents  a  distance  from  the  wall  through 
which  the  potential  outer  flow  is  displaced  due  to  the  presence  of  the 
boundary  layer.  In  other  words,  it  represents  a  mass  flow  deficit. 

Thus,  if  the  displacement  thickness  is  known  to  undergo  fluctuations  in 
time,  then  one  can  expect  mass  flow  fluctuations  near  the  surface.  These 
fluctuations  then  act  on  the  potential  outer  flow  (acoustic  medium)  in 
much  the  same  way  as  a  vibrating  surface  acts  on  its  surrounding  medium. 

The  normal  velocity  of  the  surface,  which  can  be  termed  the  piston  velocity, 
defines  a  boundary  condition  with  which  the  wave  equation  may  be  solved 
for  the  sound  pressure.  By  constructing  a  model  for  the  displacement 
thickness  fluctuations  that  occur  within  the  transition  zone,  an  equivalent 
piston  velocity  is  similarly  defined.  The  wave  equation  is  then  solved 
by  straightforward  methods  such  as  those  used  in  the  study  of  sound  radi¬ 
ation  from  randomly  vibrating  structures. 

The  concept  of  a  displacement  effect  is  not  new.  It  is  discussed  in 
some  detail  by  Laufer,  Ffowcs  Williams,  and  Childress  [A]  in  regards  to 


-6- 


27  March  1980 
GCL/pjk 


the  noise  generated  by  turbulent  boundary  layers.  The  idea  was  first 
proposed  by  Liepmann,  but  apparently  was  never  published  [4].  As  long 
as  the  wave  equation  forms  the  basis  for  the  modeling,  the  results 
derived  from  this  approach  should  be  equivalent  to  those  derived  from 
Lighthill's  analogy.  The  main  difference  is  that  the  unknown  is  now 
contained  within  the  boundary  condition  and  not  within  the  source  terms 
of  Lighthill's  equation. 

The  analysis  presented  in  this  paper  is  based  on  the  displacement 
effect.  The  power  spectrum  and  radiation  efficiency  are  predicted  using 
analytical  techniques  similar  to  those  of  Ref.  [1].  To  help  in  the  veri¬ 
fication  of  this  theory,  three  examples  will  be  presented  in  which 
experimental  data  exist. 

I.  ANALYSIS 

Following  Ref.  [1],  we  consider  the  viscous  flow  of  free-stream 
velocity  u^  over  a  flat  plate.  If  a  laminar  boundary  layer  begins  to 
form  at  the  origin  of  our  Cartesian  coordinate  system,  and  if  the  flow 
direction  is  in  the  x^-direction,  then  the  flow  will  ultimately  become 
nonlinearly  unstable  at  the  downstream  line  x^  =  Xq  (see  Fig.  1) .  There 
is  then  a  short  distance,  Ax,  over  which  turbulent  "bursts”  occur. 

Along  x^  =  Xq  +  Ax,  the  flow  becomes  fully  turbulent  while  for  x^  <  Xq, 
the  flow  is  assumed  completely  laminar. 

Within  the  source  region,  Xq  <.  x^  <_  Xq  +  Ax,  the  boundary  layer 
intermittently  changes  from  laminar  to  turbulent  regimes.  The  bursts  of 
turbulence  usually  grow  as  they  convect  downstream  at  some  mean  convection 
velocity,  u^.  The  locus  of  a  typical  burst  forms  a  wedge  of  apex  angle 
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2a,  where  we  will  assume  Emmons'  [5]  result  that  a  -  9.6°.  Now,  at  any 
arbitrary  location  within  the  source  region,  there  are  going  to  be 
instants  of  time  when  the  boundary-layer  flow  is  laminar  while  at  other 
instants  it  will  be  turbulent.  The  intermittency  factor,  y(x^),  is  a 
measure  of  the  percentage  of  time  that  the  boundary  layer  is  turbulent. 
Through  very  carefully  performed  experiments  on  a  flat  plate,  Dhawan  and 
Narasimha  [6]  showed  that  at  where  0  <  y(x^)  <  1,  the  mean  value  of 
boundary  layer  displacement  thickness  is  given  simply  by: 

«$*>  =  (1  -  y)  6*  +  y6*  ,  (1) 

where  subscript  L  stands  for  the  laminar  flow  value  and  T  stands  for  the 
corresponding  turbulent  flow  value.  We  will  assume  for  the  laminar  flow 
regimes  a  Blausius  velocity  profile;  thus,  from  Schlichting  [7]: 

6*  =  1.7208  —  Re  1/2  ,  (2) 

L  U0  X1 

where  V  is  the  kinematic  viscosity  of  the  fluid  and  Re  is  the  local 

X1 

value  of  length  Reynolds  number  (u^  x^/v) .  Within  the  turbulent  bursts 
we  will  assume  a  one-seventh  power  law  for  the  velocity  profile,  from 
which  [ 7 ] : 


6*  =  0.04625  x1  Rex  1/5  .  (3) 

Both  Eq.  (2)  and  Eq.  (3)  were  developed  under  the  assumption  of  a  zero 
pressure  gradient,  but  they  may  be  used  with  reasonable  confidence  in 
flows  with  mild  favorable  pressure  gradients. 
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The  temporal  variation  of  6  through  the  transition  zone  is  the  most 
important  ingredient  in  the  development  of  the  present  acoustic  theory. 

No  known  measurements  are  available;  therefore,  we  must  propose  a  mathe¬ 
matical  model  based  upon  physical  intuition.  This  was  the  approach  used 
to  model  the  fluctuating  wall  shear  stress  in  the  previous  analysis  [1]. 
Realizing  that  6*  cannot  instantaneously  change  from  a  laminar  value 
[Eq.  (2)]  to  a  turbulent  value  [Eq.  (3)],  we  must  introduce  a  time 
constant,  t  ,  which  mimics  the  time  required  for  6  to  change  from  6L  to 
6*  (rise  time)  or  visa  versa  (fall  time).  For  the  purposes  of  the  current 
analysis  we  will  assume  the  rise  time  equals  the  fall  time  and  that  during 
the  rise  and  fall  periods,  <5  obeys  an  exponential  dependence  with  time. 
Thus,  we  propose  that  at  a  point: 

6*  (t)  =  [1  -  ?  (t)  ]  6*  +  I(t)  6*  ,  (*) 

where  l(t)  is  a  random  indicator  function  as  shown  in  Fig.  2(a). 

Following  Laufer,  et.al.  [4],  the  linearized  boundary  condition  on 
the  acoustic  field  is 


v 

n 


Cm 


(5) 


where  4>  is  the  velocity  potential  of  the  sound  field.  From  Eq.  (4), 
Eq.  (5)  becomes 


(6) 
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* 


4 

where  A 6  =  6T  -  6^  and  I  =  9l/3t  which  is  depicted  in  Fig.  2(b).  The 

velocity,  is  the  equivalent  "piston  velocity"  of  the  surface  due  to 
fluctuations  in  the  boundary  layer  displacement  thickness.  The  acoustic 
problem  now  becomes  one  of  a  flexible  plane  surface  of  pre-scribed  normal 
velocity,  v^.  The  solution  for  the  sound  field  is  well  known.  Crighton  [8] 
gives  a  concise  derivation,  finding: 


P(r,t) 


P  3_ 
4irr  3t 


f 

|x  -  nh 

V 

n,  t  - 

JJ  n 

~  c 

S 


dn 


(7) 


where  p  is  the  acoustic  pressure,  p  is  the  fluid  mass  density,  r  =  |x| , 
c  is  the  sound  velocity,  and  S  is  the  area  covered  by  the  sources. 
Equation  (7)  shows  no  angular  dependence;  hence,  the  displacement  effect 
creates  monopole-type  sound  radiation.  This  equation  can  be  simplified 
somewhat.  Consider  for  the  surface  a  typical  length  scale,  W  (Fig.  1), 
and  for  the  transition  zone,  length  and  velocity  scales  Ax  and  u^, 
respectively.  Then,  Crighton  shows  that 


P(r ,t) 


P  9 
4frr  9t  J 
S 


v  (n,  t  -  r/c) 
n  — 


dv\ 


(8) 


if  and  only  if  W/c  <<  Ax/uc<  In  this  simplification,  the  integral  may  be 
interpreted  as  the  total  instantaneous  rate  of  mass  outflow  created  by 
the  surface. 

Substituting  for  v  ,  Eq.  (8)  takes  on  the  form: 


Ax 


p(r>  t) 


p  ±_ 

4irr  3t 


A6  (n1)t(n1,n3,  t  -  r/c)  dnx  dn3 


(9) 
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where  the  origin  of  our  coordinate  system  has  been  displaced  x^  units 
downstream. 


where  w  =  2nf  which  is  the  radian  frequency  and  T  is  the  delay  time 
associated  with  the  correlation  function.  We  would  like  to  simplify 
Eq.  (9)  one  step  farther  by  making  a  parallel  flow  assumption  over  the 

•k 

source  region.  Thus,  A<S  (r|^)  is  considered  independent  of  which  is 

k 

certainly  plausable  when  Ax  is  small.  The  term  A6  may  be  brought  outside 
the  integrals  and  evaluated  at  some  average  value  of  x^,  e.g.,  =  xq  + 

Ax/ 2. 

Equation  /10)  is  now  of  the  form: 


C(r,f)  * 


2 

9  *  9 

/(A 6  )  u) 

2  2 
8tt  r 


<t (n1 , n3 , t  -  r/c) 


i (nx  +  Zv  n3  +  C3,  t  -  r/c  +  T ) >eltl)T  dn1  dn3  d^3  dx. 
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where  I(t)  is  the  "ideal"  indicator  function  of  Fig.  3(b).  The  frequency 
depression  factor  is  easily  calculated: 


J  (ico)  = 


—  v 
2  0 


1  eiUit  dt  =  V, 


-t/t 


i  iwt  .. 
e  dt 


(13) 


The  constant  VQ  is  determined  from  the  condition  [9]  that  J(0)  =1.  We 
find  VQ  =  t^  \  from  which  Eq.  (13)  becomes: 


J(iu))  =  (1  -  ia)t^) 


and 


|j(ico)|2  =  [1  +  (U)ti)2] 


-1 


(14) 


The  integrals  of  Eq.  (12)  are  now  identical  to  those  treated  in 
Eq.  (23)  of  Ref.  [1].  We  will  not  repeat  the  details  here,  but  point 
out  that  <II'>  was  assumed  separable,  i.e.. 


<!!’>  =  R1(n1,  t)  R2(n1,  Z3)  , 


(15) 


where 


^(ip  C3)  d^3 


3  tana 


(16) 


and  R1(n1,C1,T)  =  Y^)  exp  (-a*|C1|)  exp  (-2N  |  T-C^/^  | )  (17) 


with 


a*  -  (1  +  83.35  z8)/Ax 
v  ..  .  -4.185  z2 


(18) 


4. 


I 
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N(z)  =  1.272  (uq/ Ax) ze 


z  =  n^/Ax 


-4.185  z 


The  distributions  for  the  intermit tency  factor  and  burst  frequency  given 
in  Eqs.  (19)  and  (20),  respectively  are  taken  from  Farabee,  et.  al.  [10]. 

The  power  spectral  density  of  the  far-field  acoustic  pressure  radiated 
per  unit  spanwise  width  of  boundary-layer  transition  becomes: 


3G(r.f)  P  (A<S  ^  u0  uc  *.  ,  , 

3x!""  =  - ~T~2 -  F  (kcAx’  a  Ax’  Vi/Ax’  Uc/u0)  ’ 

3  8it  r 


where 


9  ft 

.  (k  Ax)  F(k  Ax,  a  Ax,  u  /u  ) 
F*  _  C _ c _ C  U 

[1  +  (k  Ax)2(u  t./Ax)2] 

C  Cl 


.  ,  1  .  ,70  2n  -4.185z2  -4.185Z2 

F(kcAx,a  Ax,uc/uq)  =  (kcAx)2  |  - ^~~Zz  (1-~e - — - 

0  [..472  ^1  [<a*te)2+(kcto)Jl1/2 

Cl 


2a*  Ax 


[ (a*Ax) 2+(k  Ax)2] 
c 


*  * 

Yj2  ~  e  a  ^XZsin(kcAxz-(j))+e  a  ^  z  s in [k^AxQ-z )-<()]  f  dz,  (24) 


-1  *  /  / 
with  (f>  =  tan  (a  Ax/kcAx)  ,  and  kc  =  U)/uc 


9 
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Equation  (23)  has  been  integrated  numerically  by  Simpson's  rule.  Typical 

results  of  this  computation  are  presented  in  Figs.  4  and  5.  It  is  seen 

from  these  figures  that  the  magnitude  of  the  high-frequency  portion  of 

the  spectrum  depends  significantly  on  the  value  of  the  nondimens ional 

rise  time,  u^t^/Ax.  At  very  low  values  of  kcAx,  the  spectral  levels  rise 

at  12  dB/octave  and  are  essentially  independent  of  uct^/Ax.  For  kcAx  »  1, 

the  spectral  levels  roll  off  at  6  dB/octave.  The  frequency  of  the  peak 

in  the  spectrum  depends  upon  u^t^/Ax. 

As  noted  in  Figs.  4  and  5,  the  selected  range  for  the  nondimensional 

rise  time  is  between  0.05  and  1.35.  This  range  was  established  through 

the  following  dimensional  reasoning.  We  expect  that  t^  scales  with  the 
* 

length  A6  and  with  the  velocities  very  near  and  normal  to  the  surface; 

k 

these  velocities  scale  with  the  friction  velocity,  u^.  Thus,  t^  "  A6  /u^. 

Letting  u.  =  u.,/30  and  u  -0.7  u„,  we  conclude  that 
*  U  c  u 


Vl 

Ax 


Ax 


(25) 


k 

For  typical  values  of  A6  and  Ax  as  will  be  discussed  in  SECTION  II,  the 
above  relationship  suggests  that  0.05  <_  u^t^/Ax  <_  1.35. 

B .  Radiation  Efficiency 

2 

If  we  denote  the  right-hand  side  of  Eq.  (22)  by  <p^  >,  then  the 
power  radiated  per  unit  spanwise  width  is  given  by: 


N 

a 


<Pn2> 


2tt 

r 


pc 


dw 


d8 


r  sinO  d0 


(26) 
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substituting  for  <p.  >, 


we  find: 


N 

a 


*  3 

p(A6  )  uQ  uc 

2  it  c  Ax 


* 

F 


> 


(27) 


where 


* 

F 


r 

F  d(k  Ax) 
c 

0 


(28) 


The  numerical  value  of  F  depends  upon  u^t^/Ax,  where  typical  values  are 
tabulated  on  Figs.  4  and  5.  The  radiation  efficiency  is  defined  as  the 
ratio  of  acoustic  power  to  hydrodynamic  power  generated  by  the  source 
region.  It  was  shown  by  Lauchle  [1]  that  the  hydrodynamic  power  for 
transitional  flow  is  given  by: 


Nu  a  0.572  u„  Ax  a  (29) 

n  u 

where  a  =  t^Xq)  -  t^Xq)  which  is  the  difference  between  the  turbulent 
and  laminar  values  of  wall  shear  stress  at  the  beginning  of  intermittent 
flow.  Dividing  Eq.  (28)  by  Eq.  (29),  the  transition  zone  radiation 
efficiency  results: 


n 


tr 


0.278 


(30) 


Noting  that  A6  /Ax  is  only  weakly  dependent  on  u^  and  that  O  ~  pu^  , 
n  ~  Uq/c  which  is  to  be  expected  for  a  monopole  source. 
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Thus,  n  /O^dt  depends  only  upon  Re  ,  M,  and  u  t./Ax.  Figure  6  shows 
tr  1BL  Xq  c  i 

the  dependence  of  n  /rL,,,.  on  M  for  a  range  of  rise  times  and  with 
tr  i  dl 

Re  =  4  x  10^  and  u  =  0 . 8  u. . 
xQ  c  0 

Some  caution  should  be  exercised  when  interpreting  Fig.  6.  The 
ratio  of  radiation  efficiencies  presented  here  assumes  equal  source  areas. 

In  practice,  say  for  a  body  of  revolution  moving  through  water,  the  area 
of  transition  is  AxW  while  the  area  of  fully-developed  turbulent  flow  is 

AxW  which 

means  that  a  measured  radiated  noise  spectrum  may  certainly  contain 
significant  contributions  from  the  fully-developed  turbulent  flow. 

II.  EXAMPLES 

By  way  of  numerical  examples,  we  will  consider  those  cases  in  which 
flow  noise  was  measured  experimentally  under  anechoic  or  near-anechoic 
conditions.  It  is  pointed  out  from  the  outset  that  for  each  of  these 
examples,  the  experimentalists'  goal  was  to  measure  the  overall  flow 
noise  of  the  vehicle  or  fixture  at  different  points  in  space.  We  there¬ 
fore  select  those  measurement  points  which,  based  on  our  physical  insight, 
would  appear  to  be  most  influenced  by  the  noise  generated  by  the  transition 
zone.  It  is  further  pointed  out  that  many  of  the  hydrodynamic  parameters 
that  characterize  the  transition  zone  flow  are  not  available  for  the 
examples  discussed  below.  We  have  therefore  estimated  them  using  theoretical 
considerations.  The  most  significant  unknown  is  the  rise  time,  t^.  Because 
there  are  no  definitive  measurements  of  this  time  constant,  we  have  elected 
to  use  the  dimensional  analysis  that  led  to  Eq.  (25). 


of  the  same  order  as  the  body  surface  area,  A^.  Clearly, 
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A.  Case  1:  Flat  Plate  In  Air 

DeMetz  and  Casarella  [13]  performed  extensive  measurements  of  the 

fluctuating  pressures  that  occur  within  the  transition  region  of  a  large 

flat  plate  operating  in  an  anechoic  wind  tunnel.  Of  their  data,  we  select 

a  condition  where  the  microphone  was  located  within  the  laminar  boundary 

layer  so  that  its  response  would  be  due  only  to  the  noise  generated  by 

the  downstream  unsteady  flow.  A  spectrum  of  this  noise  is  presented  in 

Fig.  B-3  of  Ref.  [13].  Here,  Re  =  7.5  x  10^  and  the  measurement  took 

6  X° 

place  at  Re  =  4.758  x  10  .  Based  on  the  difference  of  these  two  Reynolds 
X1 

numbers,  r  can  be  estimated  to  be  1.55  m.  The  width  of  the  plate,  W,  is 
2.44  m,  which  is  multiplied  by  Eq.  (22)  in  order  to  obtain  the  radiated 
noise  due  to  the  entire  width  of  transitional  flow.  Now,  the  total  plate 
length  is  4.34  m  and  x^  =  4.26  m;  therefore  Ax  ^  76  mm.  In  other  words, 
transition  has  occurred  right  before  the  trailing  edge  and  the  data  should 
therefore  be  uncontaminated  by  any  noise  due  to  fully-developed  turbulent 
flow  [Eq.  (32)  predicts  Ax  >  76  mm]. 

Assuming  uc/uq  =  0.8,  the  comparison  of  our  prediction  with  the 
experimental  spectrum  is  shown  in  Fig.  7.  DeMetz,  et.  al.  [13]  attribute 
the  high  energy  levels  between  5  and  10  kHz  as  being  due  to  pinhole  micro¬ 
phone  cavity  resonance  and  plate  vibration.  Aside  from  this  band,  the 
data  compare  quite  favorably  with  the  present  model  using  uct^/Ax  =  0.15. 
Also  shown  on  Fig.  7  is  a  prediction  for  the  shear  stress  contribution 
to  the  transition  zone  noise  [1].  Clearly,  the  fluctuating  displacement 
thickness  mechanism  appears  to  dominate. 
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B.  Case  2:  Buoyant  Body  In  Water 

Haddle  and  Skudrzyk  [2]  present  flow  noise  spectra  measured  both  on 
and  far  away  from  a  buoyantly-propelled  axisymmetric  body  with  a  hemi¬ 
spherical  nose.  Because  the  terminal  velocity  is  quite  high  (30-37  knots), 
it  is  assumed  that  transition  occurs  on  the  0.48  m  diameter  body  before 
the  laminar  separation  point.  We  assume  transition  occurs  at  the  minimum 
pressure  point  which  is  predicted  to  be  at  =  0.32  m  [14].  Thus,  for 

the  37  knot  condition,  Re  =  6.27  x  10^  which  assumes  21°C  fresh  water. 

X0 

The  length  W  is  calculated  from  the  circumference  of  the  body  at  x^. 

Hydrophone  no.  7  which  is  located  at  a  45°  arc  from  the  stagnation  point 

is  selected  for  the  first  comparison.  The  estimated  point  of  transition 

is  0.13  m  downstream  from  this  hydrophone  (thus,  r  =  0.13  m) . 

The  spectrum  measured  by  hydrophone  no.  7  is  presented  in  Fig.  17 

of  Ref.  [2].  We  repeat  it  here  in  Fig.  8  along  with  the  transition  zone 

noise  predictions  for  both  mechanisms.  Again  rather  good  agreement  is 

achieved  for  the  noise  mechanism  addressed  in  this  paper.  It  is  noted, 

however,  that  the  experimental  spectrum  rolls  off  slightly  faster  than 

the  predictions.  This  may  be  due  to  a  water-borne  diffraction  loss  which 

was  discussed  in  some  detail  by  Lauchle  [3]. 

Haddle  and  Skudrzyk  also  measured  the  far-field  noise  of  the  buoyant 

unit  with  a  distant  hydrophone  array,  e.g.,  Fig.  19  of  Ref.  [2].  The 

reference  distance  is  the  radius  of  a  sphere  whose  area  is  the  same  as 

the  area  of  the  vehicle,  i.e.,  r  =  0.73  m.  This  run  was  performed  at 

30  knots,  hence,  Re  =  4.97  x  10^.  A  comparison  of  our  predictions  with 
X0 

the  experimental  spectrum  is  shown  in  Fig.  9.  The  shear  stress  contribution 
[1]  is  not  shown  because  its  directivity  function  would  predict  a  null  at 


A 
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the  beam  aspect  angle.  It  is  of  interest  to  note  that  A6  is  calculated 
to  be  0.43  mm  for  these  experiments  and  Ax  was  estimated  to  be  61  mm. 

Equation  (25)  then  suggests  that  u^^/Ax  =  0.14  which  agrees  to  within  a 
factor  of  three  of  that  value  used  for  the  prediction  which  fits  the  data 
the  best. 

III.  CONCLUSIONS 

In  this  paper  we  have  analyzed  the  role  of  a  fluctuating  displacement 
thickness  on  the  noise  generated  by  incompressible  boundary-layer  transition. 
In  essence,  the  analytical  methods  parallel  those  of  an  earlier  analysis 
in  which  the  role  of  a  fluctuating  shear  stress  was  addressed  [1].  Based 
on  numerical  calculations  of  the  power  spectral  density  and  on  the  relative 
magnitudes  of  the  radiation  efficiency,  the  mechanism  treated  in  the  present 
paper  is  found  to  be  much  more  productive  as  a  noise  source  than  is  the 
shear  stress  mechanism.  The  fluctuating  displacement  thickness  effect 
gives  rise  to  a  monopole-type  of  sound  radiation.  The  spectrum  of  the 
radiated  noise  was  found  to  rise  at  12  dB/octave,  level  off,  and  then  fall 
off  at  6  dB/octave. 

The  frequency  at  which  the  spectrum  peaks  is  not  easily  defined 
because  it  depends  upon  a  characteristic  time,  t^,  associated  with  the 
growth  of  the  turbulent  boundary  layer  in  a  turbulent  burst.  This  time 
models  the  time  required  for  the  laminar  value  of  the  boundary  layer 
displacement  thickness  to  change  to  a  turbulent  value  and  visa  versa. 

The  way  in  which  we  incorporated  t^  into  the  noise  analysis  is  heuristic. 

We  assumed  that  the  displacement  thickness  grows  from  a  laminar  value  to 
a  turbulent  value  according  to  an  exponential  time  variation  and  that  the 
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process  Is  symmetrical.  There  is  no  experimental  evidence  that  this  is 
true  (or  for  that  matter,  any  other  functional  form.)  Based  on  the 
comparisons  of  this  theory  with  experimental  data,  particularly  those  of 
Case  1,  it  appears  that  these  assumptions  may  be  realistic. 

A  basic  problem  to  be  resolved  deals  with  the  measurement  of  t^. 
According  to  Laufer,  et.  al.  [4],  <S  (t)  can,  in  principle,  be  measured 
in  a  laboratory  frame  of  reference  by  performing  simultaneous  measurements 
of  the  fluctuating  velocity  and  pressure  and  their  correlation  through 
the  boundary  layer.  One  then  couples  these  data  with  a  first  order  partial 

•k 

differential  equation  for  6  derived  by  integrating  the  non-steady 
continuity  equation  across  the  boundary  layer.  If  this  procedure  can 
be  demonstrated,  then  it  is  conceivable  that  t^  can  also  be  measured. 

Such  measurements,  to  this  author's  knowledge,  have  never  been  carried 
out.  Perhaps  the  findings  presented  in  this  paper  may  provide  the  impetus 
required  for  future  experimental  investigations  of  this  type. 
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Schematic  representation  of  the  transition  from  laminar-to- 
turbulent  flow  on  a  flat  plate. 
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Fig.  2.  The  indicator  function  (a)  and  its  first  time  derivative  (b)  which 
is  used  to  model  the  random  variation  in  displacement  thickness 


during  transition. 
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Comparison  of  the  theory  with  an  experimental  spectrum  measured 
on  a  flat  plate  by  DeMetz  and  Casarella  [13]. 
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Fig. 


.  Comparison  of  the  theory  with  an  experimental  spectrum  measured 
on  the  surface  of  a  buoyantly  -  propelled  vehicle  by  Haddle  and 
Skudrzyk  [2]. 
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